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We study the destabilization of a shear layer, produced by differential rotation of a rotating 
axisymmetric container. For small forcing, this produces a shear layer, which has been studied 
by Stewartson and is almost invariant along the rotation axis. When the forcing increases, in- 
stabilities develop. To study the asymptotic regime (very low Ekman number E), we develop a 
quasi-geostrophic two-dimensional model, whose main original feature is to handle the mass conser- 
vation correctly, resulting in a divergent two-dimensional flow, and valid for any container provided 
that the top and bottom have finite slopes. We use it to derive scalings and asymptotic laws by a 
simple linear theory, extending the previous analyses to large slopes (as in a sphere), for which we 
find different scaling laws. For a flat container, the critical Rossby number for the onset of insta- 
bility evolves as E 3 ^ and may be understood as a Kelvin-Helmoltz shear instability. For a sloping 
container, the instability is a Rossby wave with a critical Rossby number proportional to fiE 1 ^ 2 , 
where p is related to the slope. We also investigate the asymmetry between positive and negative 
differential rotation and propose corrections for finite Ekman and Rossby numbers. Implemented 
in a numerical code, our model allows us to study the onset over a broad range of parameters, 
determining the threshold but also other features such as the spatial structure. We also present a 
few experimental results, validating our model and showing its limits. 



I. INTRODUCTION 
A. Background 

The destabilization of a shear layer in a rotating system is of general geophysical interest. Shear layers dominated 
by the Coriolis force, are found in natural systems, as in the atmosphere or in the ocean, as well as on other planets 
like the jets of Jupiter, which might be unstable and responsible for the observed big eddies. 

This work considers the shear layer produced by differential rotation of the boundaries of a fast rotating container, 
in axisymmetric geometries. The linear problem resulting in cylindrical shear layers aligned with the axis of rotation 
has been studied by [2j| in the spherical shell geometry. He pointed out the nested shear layer structure and their 
scaling. This problem is closely related to the simpler case of a split cylinder solved earlier by too. In that case, 
the structure of the shear layer is simpler, but there is still a nested structure. The general structure of these detached 
shear layers is studied bv |l7j by focusing on the singularities of the bound ary-layer equ ations. Later, a numerical 
study has been done by Q, recovering the features and scalings predicted bv IStewartsonl 

The linear stability of this kind of layers has been theoretically investigated by Q who applied a linear stability 
theory to the pressure equation in both cases of constant and varying depth. He neglected the bulk viscosity but kept 
the Ekman friction, recovering some kind of Rayleigh criterion, showing that for the flat case, one should expect a 
critical Rossby number Ro c ~ _E 3 / 4 (where E is the Ekman number). For a varying depth container, he discussed the 
criterion obtained by 0| . 

One may notice that the difference between flat and varying depth containers has been made explicit by Q for 
thermal instabilities in rapidly rotating systems : for a flat container the instability is of Bcnard type (critical Rayleigh 
number Ra c independent of Ekman number), whereas for a spherical shell the instability is a thermal Rossby wave 
(Ra c ~ £?-V3). 

The setup of [HI, with a central horizontal disk rotating differentially in a flat cylindrical tank, gives unexpected 
results which are not yet well understood. They also observed a strong asymmetry between positive and negative 
differential rotation. Later, [Tsj performed a similar experiment but with the differential rotation concentrated on the 
bottom of their flat tank. They find a critical local Reynolds number for the destabilization of the flow, independent 
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FIG. 1: Critical Rossby number versus Ekman number for the spherical shell geometry. The results obtained numerically 
by H3 and our lab experiment are in good agreement. The symbols represent established instability, and the error bars the 
domain where instability was not always clear. "3D simul" is a run with our own 3D code. 



of the Ekman number. Their linear theory which includes both Ekman friction and bulk viscosity is in very good 
agreement with their results. They point out that the bulk viscosity is very important to reproduce accurately the 
experimental results. They showed that a critical local Reynolds number (defined for the shear layer) describes the 
threshold of instability properly. 

Later, reproduced almost exactly the case studi ed bv |2^|. The s tability threshold they found is still governed 
by a critical local Reynolds number, but the theory of lNiino fc Misawal predicts a value much lower than the observed 
one. They also observed that the behavior of the flow does not depend strongly on the sign of the differential rotation. 

The shear- flow instabilities in a parabolic tank and a free surface was investigated by |lal[ . They studied both cases 
with constant depth or varying depth and pointed out that the varying depth has a stabilizing effect, as predicted by 
the linear theory. However, their setup uses the centrifugal force to adjust the slope of the free surface, so that the 
slope is directly related to the rotation rate. 

More recently, using a 3D code, ^3 studied the differences between positive and negative differential rotation for 
the onset of these instabilities, pointing out the importance of the geometry of the container. 

[T^ l computed numerical solutions for the instability of the Stewartson layer between two rotating spheres. Mean- 
while, we ran the corresponding experiment, measuring the threshold. Details of the experimental setup can be found 
in WI C II The results of theses two studies are in agreement, and summarized in figure [T] The study of showed 
that the asymptotic regime of the Stewartson problem cannot be illustrated numerically for Ekman numbers higher 
than 10 -5 . Unfortunately, neither the experimental nor the numerical approaches were able to reach Ekman numbers 
lower than 10 -5,5 . This number should be compared with geophysical flows, where Ekman numbers can be much 
smaller : for the Jovian atmosphere, the oceans and the Earth core the Ekman number is smaller than 10~ 10 . 



B. Outline 



To study the very low Ekman number regimes we use a quasi-geostrophic (QG) model, detailed in iflTl : we derive the 
equations of an improved QG-model, that will be used to address the problem in the next sections, through numerical 
computations and theoretical analysis of the equations. This QG-model is improved over previous ones by taking into 
account the mass conservation and Ekman friction. The latter is also dynamically important for non-linear regimes 
because of its dissipativc nature. We were able to benchmark this model with some fully 3D computations by 01 
and the results are quite encouraging. 

flllll presents results concerning the Stewartson layer, comparing the QG-model with existing results, showing its 
ability to reproduce the axisymmetric flow. A numerical code based on the quasi-geostrophic approximation allows us 
to reach Ekman numbers as low as 10 -10 which is very likely to lie in the asymptotic regime, from which extrapolations 
can be made to even lower Ekman numbers. 

In fllVI using numerical calculations as well as theoretical arguments, we investigate the spatial structure of the 
instability and its asymptotic evolution when lowering the Ekman number. When a /3-effect is present, the instability is 
a Rossby wave, and consequently we find that its spatial structure changes dramatically with the sign of the differential 
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rotation. This explains the asymmetric sensitivity to geometry studied by [12j. Another interesting feature is that 
the radial extent of the instability may be independent of the Ekman number in some particular cases. 

In f]^we derive the scaling laws expected in the asymptotic regime, focusing on the stability threshold. For a flat 
container, an i? 3 / 4 law is recovered for the critical Rossby number. We extend this study to finite slopes, in which 
case the scaling changes to /3E 1 / 2 , where (3 is related to the slope. This scaling is the one that should be applied to 
spherical shells. We perform an analysis leading to a Rayleigh-like criterion, which allows us to evaluate quantitatively 
the effect of the asymmetry induced by the Rossby wave mechanism. Numerical results obtained with our QG-model 
are also reported here, for various geometries. 

WII compares our results with previous experiments, as well as with our own, to span a wide range of different 
geometries. Despite the lack of quantitative agreement between numerical results obtained with the QG-model and the 
experimental data, the global behaviour of the instability is captured. Finally we conclude in Will with a discussion. 



II. QUASI-GEOSTROPHIC MODEL 

The two-dimensional quasi-gcostrophic (QG) model we present here has been developed to study the dynamics of 
a fluid in a rapidly rotating container, invariant by rotation around the rotation axis, without being restricted to 
the small slope case. We use the cylindrical coordinate system (e r ,e,p,e z ), with r the distance to the axis, </> the 
angle in the plane perpendicular to that axis, and z the height. The container is defined by two functions z = £+(r) 
and z = L_(r), for the top and the bottom surfaces respectively. Here, we first restrict to the simpler case where 
L{r) = L + {r) = -L r {r). 

In a frame rotating at angular velocity Qq around the vertical axis z, the flow is described by the Navier-Stokes 
equation for an incompressible fluid with constant density, including the Coriolis force : 

du 

— + (it. V) u + 2O A u = - Vn + vV 2 u (1) 

where u is the fluid velocity field, v the kinematic viscosity, and II the reduced pressure field, including centrifugal 
and gravity potentials. We introduce a velocity scale U, a time scale Slo - 1 and as typical length scale R the radial 
size of the container. We then rewrite equation ^ with non-dimensional quantities : 

— + Ro (tt.V) u + 2e z A u = -VII + E\7 2 u (2) 
dt w 

with Ro = U/RHq the Rossby number, E = v/R 2 Q,q the Ekman number, and e z the unit vector along the rotation 
axis. 

In our study, we take for typical velocity the one based on the differential rotation U = RAD,, which is supposed to 
dominate. Hence, the Rossby number reduces to Af2/f2 . 



A. z-averaged vorticity equation 

For vanishing Ekman and Rossby number, the Proudman- Taylor theorem ensures that the flow is invariant along 
the rotation axis direction = 0). This suggests that for small Rossby and Ekman numbers, the z- variations 
should be weak and we expect the flow to be mainly two-dimensional for long-period motions (3> fig" 1 )- 

We define the z-averaging ( ) by 

1 f +L 

(u) = — J u(r, z) dz (3) 

Keeping in mind that L depends on r, we have 

d (u) I du\ 1 dL ( u(r 7 +L) + u(r, —L) . , .. 

{u{r,z)} 



dr \dr / L dr 



Hence, a sufficient condition for the r-derivative to commute with z-averaging is that u is independent of z, that is 
u = (u). 

We now average the z-componcnt of the curl of equation [2] with the additional constraint that u r and are 
independent of z (so that we can permute r-derivation and z-averaging). The z-component of the vorticity u> is 
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denoted by u and is also z-invariant. We obtain 



Note that we keep the vortex stretching term (^f-)- All other terms vanish, without any condition on it 2 . 

B. Mass conservation 

The local mass conservation equation for an incompressible fluid, 

1 dru r 1 duri, du. 



r dr r dtp dz 



(5) 



implies that is independent of z, so that we can now draw a global picture of the flow in the frame of our 
QG-model. 

• u r and Ufj, are z-invariant 

• u z depends linearly on z. 

C. Boundary conditions and Ekman pumping 

At the boundaries, Ekman layers will smooth out the vorticity jump. Their thickness is supposed to be small 
(E <C 1), and we will use the asymptotic expression given by |lOj |. equation 2.6.13 : 

u\n\ ±L = lE^n.V A ( "A(i»-yT(n-^ (6) 



V 7 ! 



n.e. 



where Ub is the velocity of the boundary. An important point is that expression is valid for time-dependent flows 
with time scales larger than a few rotation periods. Greenspan used that expression to study the spin-up process in a 
sphere, which is a time-dependent flow with time scale of order E 1-1 / 2 , and Q show experimental data in very good 
agreement with this linear theory. It is the suitable approximation for our study, as we will find typical period of 
order E^ 1 / 4 in section This gives us the appropriate boundary conditions for our problem, allowing us to express 
the vertical velocity. Defining (3 



L dr 



(7) 

z=L 



we have 

du z 
dz 



E 1 / 2 .P(u r ,u^,r)+f3u r (8) 



where P is a function depending on the geometry of the container and linear with respect to the jump of u r and at 
the boundary and to their first order derivatives. See Appcndix^for the explicit expression for a spherical container. 
In the case of a constant-depth container, equation [S] reduces to 

dv E 1 / 2 
dz 2L 

where Of, is the vorticity at the boundary. 

In the case of vanishing Ekman number, the Ekman pumping could be neglected compared to the (3 term (see AC 2|) . 
However, the Ekman friction is a dissipative term whereas the /3-term is not, it is thus expected to play an important 
role for the non-linear dynamics of quasi-gcostrophic flows. 

The problem is now reduced to the determination of two components of the velocity, u r and u^, depending only on 
r and tp. We also need to set a boundary condition on the sides of the container, to be consistent with the modeling 
of (^f-)j we have to use a no-slip boundary condition u r = = 0. At the center (r = 0) we impose that u r and 
remain finite. 
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D. Scalar pseudo-stream function 

We introduce a scalar field tp, denned by 

I dip 

* S (10) 
Uff, is now constrained by the three-dimensional mass conservation equation^ so that for non-axisymmetric flows, 



» = -W-0{-£ + 7)-w* (12) 

With these equations, we can see that the /3-terms are negligible for /3A <C 1, where A is a typical radial length scale of 
the velocity field variations. The model then reduces to the widely used small slope approximation [see e.g.HLI3.ll^. 
It will be shown later, that at the stability threshold, the typical length scale is of order E 1 / 4 , so that asymptotically 
and at the threshold of instability, the small slope approximation may even be valid in a sphere. However, increasing 
the forcing, large scale structures will eventually appear, with up to order-one scales, so that we cannot neglect the 
/3-terms any more in a fully turbulent computation. 

Similar approaches have been developed recently. |l4j | proposed a formulation for ocean modeling but without 
the mass conservation correction. include the effect of weak topography and Ekman friction. The advantage of 
their formulation is that the Ekman pumping is considered as a source of horizontal divergence. However the simple 
expression they used for the Ekman boundary condition implicitly assumes nearly flat topography (small (3). 

E. Axisymmetric flow 

For the axisymmetric flow, the scalar function formalism cannot be applied. As it is simpler to work with the velocity 
field rather than with the vorticity, we project the Navier-Stokes equation [2 on e^, noting that (e z A u) .e^ = u r and 
remembering that = 0. Averaging over <j> to keep only the axisymmetric component, we obtain 



^ + «o( U ^ + ^ +2W,. = £ (v%-lf) 



where u stands for the </>- aver age of u. 

The mass conservation equation reduces to 



1 dr u r du z 
r dr dz 







and using equation |H1 we can show that u r is only due to the Ekman pumping and is directly related to u$ so that 
the axisymmetric part is actually one-dimensional : u = u^ir) + 0(E 1 ^ 2 ). However, non-linear interactions of 
non-axisymmetric modes can produce axisymmetric flow, so that we keep the terms involving u r . Furthermore, the 
Ekman pumping flow u\ of order E 1 / 2 cannot be neglected. Equation reduces to 



with 



i x r =E^—{u^-U b ) (14) 



(15) 




F. Numerical implementation 



We wrote a numerical code, that solves the equation of the model described above using the pseudo-stream function 
formulation. We use a finite difference scheme in the radial direction, and a Fourier expansion in azimuth. The time 
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FIG. 2: Sketch of the nested Stewartson layers. The arrows show the Ekman pumping inside these layers. 



integration is done using the semi-implicit Cranck Nicholson scheme for the linear part, and a simple Adams-Bashforth 
scheme for the non-linear terms, which are computed in direct space, as in 

Our numerical implementation is able to handle various geometries from flat containers to spherical ones. 

When compared to fully 3D computations, our QG-model gives quite good results for the basic state flow (see SHUf) 
as well as for the stability threshold (see fig. 1121 and WIBI for more details). 



III. STEWARTSON LAYERS 



The linear studyof shear layers due to differential rotation of the walls has been done in the case of two concentric 
co-rotating disks j^, and in the case of two concentric differentially rotating spheres |23j . 

At first order, the flow is invariant along the direction of the rotation axis, and consists of a shear layer which is 
smoothing out the velocity discontinuity. There are actually two nested layers (see fig. EJ with width scaling like E 1 ^ 3 
for the inner layer, and E 1 / 4 for the outer one. Note that the 1/3 layer is not two-dimensional and smoothes out the 
singularity in the Ekman pumping due to the outer layer. For the spherical shell case, the outer layer width scaling 
is slightly modified to E 2 ! 1 for r < ro only. Thus, we can expect that the destabilization of the shear layer in both 
cases will be closely related. 



A. QG-model predictions 

In the linear and stationary regime, when the Rossby number is small equation 1131 reduces to a balance between 
the Ekman pumping term ll4l and the viscous forces in the bulk. 

2El/2 ^(^-U b ) = E(v 2 u^^ (16) 

which defines the Stewartson l ayer. Introducing A, the width of the Stewartson shear layer we easily recover the 
scaling obtained bv lStewartsonl 

A ~ E 1 / 4 (17) 

Note that the neglected non-linear terms are small if Ro/A < 1. In JVAl we will find that Ro <~ E a with a > 1/2, 
validating this approximation near the threshold. 

One of the major issues with the QG-model, is that it cannot model the inner E^Mayer. However, according to 
the numerical study of Q, the £' 1 / 3 -layer has a negligible effect on the velocity profile at Ekman numbers smaller 
than E = 1CT 6 . 



B. Numerical results 



The split-sphere geometry consists of a sphere, split at a given cylindrical radius r$ (see fig. |3J). We solve equation 
1161 in the frame rotating with the outer part at angular velocity while the inner part rotates differentially at 
Af2 (positive or negative). This linear problem is solved by a simple matrix inversion using a finite difference radial 
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FIG. 3: (a) the spherical shell geometry, studied by p3|. (b) the split-sphere geometry this paper focuses on. 




FIG. 4: Basic state flow showing the shear layer which scales with E 1 ^ 4 in the asymptotic low Ekman number regime, for 
ro = 0.35 in the split-sphere geometry, (a) shows the raw radial profiles and (b) the same profiles rescaled radially, which merge 
together when E is sufficiently low (E < I0~ 6 ). 

scheme with variable grid spacing. In this case, the boundary velocity Ub is a simple function, Ub = r if r < ro and 
Ub = if r > ro. The resulting shear layer for ro = 0.35 is shown on figure^] Properly scaled, all these profiles can be 
superposed in the asymptotic, small Ekman number, regime. The size of the shear layer scales like E 1 / 4 , as predicted 
by [221 and the asymptotic regime seems to be attained for E < 10~ 6 . 

The spherical shell (see fig. |3J geometry is the one studied originally by [2(]], and is also of geophysical relevance 
because the core of the Earth is composed of a solid iron ball (ro = 0.35) surrounded by liquid iron. The structure of 
the Stewartson layer in a spherical shell (including the £' 1 / 3 -layer), has been modelled numerically by 0, with a fully 
three dimensional code. We have compared their results with the ones obtained with our QG-model for the spherical 
shell geometry for Ro = and find that the difference between the two basic state profiles of the angular velocity is 
always smaller than 0.03 for E = 10~ 6 ; and for E = 10~ 8 the difference never exceeds 0.008. Hence, we may expect 
the .E^-layer to play only a minor role at such small Ekman numbers. 

IV. SPATIAL STRUCTURE AT THE ONSET 

Using our QG-model, we obtained numerically the first unstable mode for four geometries : the flat container; 
a container with an exponential shape (ft = — 1 everywhere); the split-sphere geometry (see fig. 01; and the disk 
geometry. The latter is a split-sphere where the polar caps have been replaced by flat disks, to mimic our experimental 
setup (see fig. IT3|) . The first unstable mode is represented on figure[S]for all these geometries. Note that in the reference 
frame chosen for our computations and for small enough Ekman numbers, the mean flow vanishes in the outer region, 
whereas the inner region is almost in solid body rotation at the angular velocity Af2 = Rotto (see figure 0J. 

A. Dependence on the sign of the Rossby number for j3 7^ 

We can observe on figure that many features depend on the sign of the Rossby number, for /3 ^ 0. The most 
striking one is the spatial behaviour of the split-sphere arrangement (fig. ISJd) : for Ro > the instability fills the 
outer region (r > ro) with spiralling arms, whereas for Ro < the critical mode seems to extend more in the inner 
region (r < ro), leaving the outer region almost at rest. The constant (3 case (fig. |5]a,) also shows long range extension 
in the outer region for Ro > 0, but no spiralling. 

To understand all these features, we must turn to the instability mechanism, which is strongly dominated by the 
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FIG. 5: North-pole views of the radial velocity for the first unstable mode at E = 10 for various geometries. Red is positive, 
blue negative, and green zero. The picture on the left is for Ro > 0, and the one on the right is for Ro < 0. On each picture, 
(a) shows the constant /3 = —1 geometry; (b) shows the split-sphere case; (c) is the disk setup; and (d) is the unstable mode 
for a flat container (,3 = 0). 

U Shallow , '^veiling wave 

L> Ro>0 deep 

U shallow travelling wave 

Ro<0 deep n 

FIG. 6: Instability mechanism showing a different spatial behavior depending on the sign of the Rossby number 

slope p. FigureElsketclies this mechanism. For a column of local vorticity lo moved to a deeper region, the conservation 
of mass and total vorticity (2f2o + ui > 0) leads to the following behaviour : If w > the local vorticity is amplified, 
and if w < it is damped. 

Hence, for Ro > the vorticity sheet is negative implying that the perturbation is damped in the deeper region, 
whereas it is amplified in the shallow region. For Ro < 0, it is the opposite. In both cases, the amplified vorticity will 
induce a flow which induces a drift of the pattern in the prograde direction (as for any Rossby wave) . 

This fundamental asymmetry between positive and negative Rossby numbers partly explains the results obtained 
by We will discuss an example in WI Bl 

B. Radial scalings 

In figure [3Jtbc the perturbation of the shear layer spreads across the whole volume whereas in the flat container 
case (fig. |SJi) or the thermal convection instability in a sphere [see@, it is radially localized. This is explained by a 
simple model detailed in appendix [51 and the main results obtained are summarized in the following lines. 

For small (3 as in the inner part of the split-sphere arrangement, the critical mode has a radial extension comparable 
to its azimuthal scale, without oscillations; but for large (3 as it is the case for the outer part of the split-sphere 
arrangement, the instability consists of a rapid oscillation of radial length scale A related to the azimuthal wave- 
number to by 



A - myJ\2P\ 



(18) 



0.5 0.55 0.6 0.65 0.7 0.75 0.8 0.85 0.9 0.95 1 

FIG. 7: Vorticity profile at fixed <j> at the stability threshold for Ro > in the outer region. 




FIG. 8: Evolution of the radial size (evaluated using spatial derivation) of the unstable mode near the equatorial singularity for 
the split-sphere geometry. All these curves exhibit a wide range of slope 1/2, which is the expected trend to cancel phase-mixing. 
Near the equatorial boundary (r = 1) this scaling is no longer valid. 



This rapid oscillation is contained in an envelope of radial size 

A- (2\/3\y 1/2 a' 1 E° (19) 

Eg uat ion 1191 implies that the radial extent of the perturbation is independent of E. This spatial structure corresponds 
to the spiralling of the Rossby wave observed for variable (3 (see fig. Etc). We may also notice a fundamental difference 
with the thermal convection case investigated by |24| and jl3l | where the E 1 / 6 radial extension decreases with E, so 
that the instability remains localized for small enough E. 



1. Numerical results 



We show in figure [7| a radial vorticity profile (at a fixed 4>) obtained with the split-sphere setup. It shows that 
the vorticity gets amplified by the decreasing depth as approaching the equator, while the spatial frequency of the 
oscillations increases, in agreement with equation^] In addition even when decreasing E, the perturbation seems to 
spread across the whole width of the tank, as predicted by the scaling ^| This is not due to the singularity of the 
QG-model at the equator, because we also observed the same phenomena with a truncated sphere, where the slope 
remains finite (see flV C 111 . 

The scaling 1181 is observed quantitatively in numerical simulations as shown on figure |H1 On this curves we can 
also see the effect of the equatorial boundary, which modifies the scaling near r = 1. When (3 is constant, there is no 
spiralling but the instability still fills the whole gap, as figure shows. 
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V. ASYMPTOTIC LAWS FOR THE ONSET 
A. Scaling analysis 

We perform a linear stability analysis on equation 0] Let U(r) be the basic state velocity profile solution of the 
linear equation [TBI Cfor Ro = 0), and u the perturbation. With the vorticity associated with U(r) : 

_ dU U 
tt = — + - 
or r 

we find 

_ / dn U duj\ „ / du z \ . . 

W + i?0 (^ + "^)" 2 \^/^ V " (20) 

where we have neglected the small local vorticity Ro Q. We now introduce the frequency / of the disturbance, and 
its azimuthal wave-number m. 

At the stability threshold, we expect the inertial terms to be both of the same magnitude, 

A" 2 ~ m 2 

which implies the wave-number of the perturbation to be comparable to the thickness of the shear layer 

1. Without slope 

When there is no slope equation |S] reduces to 

'du z \ E^ 2 
so that equation 1201 gives in terms of order of magnitude : 

E 1 / 2 

mf + Ro A~ 2 + Ro m 2 — m + Em 3 

L 

With m ~ E- 1 ^, the stability threshold is given by a balance between the non-linear forcing and the viscous damping 

Ro c ~ E 3 ' A 

which is the result obtained by the linear theory of Q and by the critical local Reynolds number (Re ~ RoE^ 3 / 4 ) 
theory of ■ ^ is a ls° consistent with the experimental results of an d an d corresponds to a Kelvin-Helmoltz 
instability of the shear layer. 

For an oscillatory instability, the frequency / should be consistent with mf ~ Ro m? , and so 

This effect, if present, may be very difficult to separate from the advection of the pattern, because they are of similar 
amplitudes. 

2. With slope 

When the slope is important, ^f- ~ /3u r , and so equation 1201 gives the following balance : 

mf + Ro A- 2 + Rom 2 + 2/3 ~ Em 3 
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FIG. 9: Basic state axisymmetric flow profiles (solution of eq. 116^ . as a function of radius r for E = 10 6 and a split-sphere 
geometry with ro = 0.35. (a) rotation velocity U/r; (b) vorticity Q; (c) radial derivative of vorticity 



If again m ~ E L ' 4 , the viscosity is negligible and the stabilization is due to the slope : the instability has to overcome 
the Proudman- Taylor constraint. As a matter of fact, from a balance between slope and forcing, we find 

Ro c ~ 2(3E 1/2 (22) 

This is the expected scaling for non-flat containers such as spherical shells, and is quite different from the flat case. 
This result agrees with the criterion obtained by 0] applied to a Stewartson layer. 
For a time-dependent instability, the frequency / verifies mf ~ 2/3, and 

which is the Rossby wave dispersion relation. At the onset, we may notice that the Ekman pumping is comparable 
to viscous dissipation in the bulk of the fluid. 



B. Rossby wave asymmetry effect 



As stated in flIV Al the instability develops as a Rossby wave. Here we give a Raylcigh-like criterion for the onset 
of instability and use it to evaluate the asymmetry induced by the Rossby wave. 

The basic state flow U, represented on figure is the solution of the linear equation 1161 The sign and amplitude 
of the forcing are both controlled by the parameter Ro in equation [201 We are looking for solutions of the form 
ip(r,(f),t) = e ml ^ +xt ip(r). We use the tp formulation introduced in section Hi Dl but without the /3-terms to keep 
simpler expressions. Ignoring Ekman friction and viscosity, we obtain from equation I20I 



2/3 - Ro^- ) 



Multiplying by ip* L(r)r and integrating over r we have 



lo [ip] ip* Lr dr 



1 2/3 



or 



Ro— 



im L dr 



X + mL RoU 

Expanding the left hand side of this equation, and integrating by part, we obtain 



1 d 



dip 



1 m 2 



i>* —^i> r— Lr dr = / -—ipip*Lr dr + 

\ r z r or \ or J J Jo r 

and from the boundary conditions the bracketed term vanishes, leading to 



Or 


l 


/' 


dip 





Jo 


dr 



Lr dr 




dip 



dr 



Lr dr 



.M 2 (2ff- Ro^) (\*~ ifRoU) 
\\ + «zlr U\ 2 



Ldr 



With 5R(A) the real part of the complex number A, the imaginary part of this equation gives 

nl (2/3 



= 5R(A* 



Ldr 



\X+i^RoU\ 



(23) 



(24) 



(25) 



(26) 



(27) 
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For a mode to be unstable, we need 3?(A*) 7^ 0, which means that there is an r for which 

dfl 

2(3 - Ro — = (28) 
or 

This is exactly the criterion already obtained by |l6[ . Using the complete pseudo-stream function given in section 
III Dl would lead to the same result. Therefore, the determination of the stability threshold will not be affected by 
using the small slope approximation, as it is usually done for the thermal convection case [seeQ. 

The critical Rossby number Ro c that will satisfy [2*51 will be obtained for the maximum value M of |^r|- We must 
now take care of what value of (3 comes into that expression. Actually, the shear layer samples the values of (3 about 
A to the left or to the right of the split radius ro, depending on the sign of Ro. Then we have 

11 M y ' 

At that point, remembering that the £^ 1 / 3 -layer smoothes out discontinuities in ^p, we can expect that it will reduce 
the value of M, resulting in an higher Ro c . Hence, the effect of the £' 1 / 3 -layer may be essentially to increase the 
critical Rossby number. 

Assuming (3 is a smooth function and A is small enough to use a Taylor expansion 



|/J(r ±A)| = |/3 |±A 

we obtain 

\Ro ± \ =2 A 2 (\(3 Q \ ± A 



d\ 



Or 



This shows a fundamental asymmetry between positive and negative Ro due to the Rossby-wave nature of the 
instability. 



C. Numerical results 

1. Split-sphere geometry 

In the split-sphere geometry, the slope parameter (3 becomes very large near the equator and the expression of the 
Ekman pumping we use is no longer valid near the equator as the Ekman layer becomes singular f |2l|: see also flA 3|l . 
However, this area is very small and so we expect that it will not affect the determination of the stability threshold. 
To test this we followed the idea of [25| and truncated our container at r = 0.9 to prevent the slope from growing 
to large near the equator. The solutions showed the same features as the untruncated calculation and very similar 
stability thresholds. 

The results for the stability threshold are given in appendix IC II and plotted in figure ITUl They seem to match the 
asymptotical scaling obtained previously, especially for small Ekman numbers. We expect the perturbation to be a 
Rossby wave and in order to compare its pulsation with the asymptotic laws, we need to correct for the advection by 
the mean flow. Considering that the outer region (where the Ro > instability develops) is not advected, we don't 
correct the Ro > data. However the inner region (where the Ro < case develops) can be considered as in solid 
body rotation (for small Ekman numbers) and the simplest way to correct for the advection is then to subtract the 
rotation rate Af2 from the phase speed c = u>/m for Ro < 0. The pulsations of the Rossby waves are plotted in figure 
I1UI 1. and they agree quite well with the asymptotic law, for both positive and negative Ro. 

We checked that the difference between critical positive and negative Rossby number is proportional to i? 1 / 4 , in 
agreement with the findings of H V HI 

With our numerical code, we are able to test various boundary conditions. We report that the boundary conditions 
(free-slip or no-slip) hardly affect the stability threshold, whereas the presence of Ekman friction has a small yet 
visible effect. The divergence correction we introduced in All Dl seems to be negligible at the threshold. The results 
of these tests are summarized in table IC 21 
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FIG. 10: Stability threshold data (given in table ICTt as a function of the Ekman number, (a) represents the threshold; (b) 
shows the critical wave number; (c) is the raw pulsation at the threshold; (d) is the pulsation corrected for advection effects. 
Dotted lines indicate the asymptotic exponent obtained in § V A 21 
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FIG. 11: Critical Rossby number (a) and critical wavenumber (b) for various geometries. Dotted lines indicate the asymptotic 
exponent obtained in W A 21 



2. Other geometries 

We report the results of the numerical computations for the four studied geometries. To compare the flat and 
non-flat case, we set up different geometries and do all the computation using our numerical implementation of the 
QG-model described in SjH] For the flat geometry we choose a flat cylindrical container with a height equal to its 
diameter. We also set up a container with an exponential shape, defined by (3 = — 1 everywhere. The split-sphere 
geometry studied in detail in W C II and finally the disk setup. This latter geometry is a split-sphere where the polar 
caps have been replaced by flat disks to mimic our experimental setup (sec fig. II 31) - However, with the QG-model 
we are not able to take the height-jump at the edge of the disks into account, so the resulting disk geometry is in 
fact flat for r < r$ and spherical for r > rrj, with no height jump. The results obtained are summarized in figure ITT1 
The critical wave-number data shows that the evolution of m c follows an E 1 / 4 law, for both flat and non-flat cases, 
as expected from the asymptotic study ( j}V All . For the flat geometry, there is no difference between positive and 
negative Rossby number. The split-sphere setup, the j3 = — 1 case and the disk geometry for Ro > have comparable 
critical wave-numbers. 
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E i?of mf i?of mf 
1.1 10~ 5 0.044 6-8 0.021 6 
8.3 10 -6 0.032 6-8 0.017 6 
5.6 10" 6 0.025 6-8 0.013 7 

TABLE I: Comparison between the experimental results of (denoted by . F ) and the numerical results we obtain (. s ) for the 
stability threshold. 

E Ro? Rof mf 

3.0 10~ 5 0.050 5-6 0.048 5 
6.0 10" 6 0.015 7-9 0.0145 8 

TABLE II: Comparison between the experimental results of [T3 | (denoted by . N ) and the numerical results we obtain (. s )for 
the stability threshold. 

The flat geometry has much lower critical wave- numbers and Rossby numbers. From figure we also see that the 
first unstable mode is totally symmetric when changing the sign of Ro. It is also different from all the other cases 
(a,b and c) which present important slopes. The disk geometry seem to be a hybrid case between flat and non-flat 
geometry : From Eh, the spatial structure for Ro > is very much like the split-sphere setup, whereas for Ro < 
it seems to be somewhere in between split-sphere and flat geometry The figure lllb confirm this, with a critical 
wave-number between flat and non-flat. This is because the Ro < instability develops in the flat inner region. On 
the other hand, figure ITTk shows a critical Rossby number for Ro < very close to the split-sphere arrangement, 
whereas the Ro > seems to follow rather a slope-like E 3 ^ 4 law, even if the threshold is significantly higher than for 
the flat container. 



VI. COMPARISON WITH EXPERIMENTS 
A. Flat container 

We reproduce numerically the constant depth experiments done by @ . The results of our linear calculation for the 
threshold are significantly lower than the experimental ones, but are in good agreement with the theory of |l8j . The 
critical wave-length is in good agreement as shown in table [I] 

To explain the discrepancy between experiment and theory, || invoke the lack of the inner E 1 ^ 3 layer in the 
modeling of which is also missing in this QG-modcl. 

We also reproduce numerically the experiments of |l8j| , which has the particularity of having differential rotation 
only at the bottom of the tank, while the top is rotating at the mean angular velocity. Results are reported in table 
[n]and agree surprisingly well with the original experimental data, as did their numerical calculations. 

B. Spherical shell geometry 

The case of a spherical shell is quite singular : (3 has opposite sign for r < ro and r > ro, and the magnitude of 
(3 is much higher for r < ro. This implies that Rossby waves travel in opposite directions in these two regions. For 
positive Ro, the instability may develop either in the r > ro region or in the r < ro one (see i flV Al and fig. [BJ . In 
both directions the depth is decreasing, but from eq. 1221 we know that the instability will appear in the outer region 
(r > ro) where the magnitude of j3 is lower. However, for negative Ro there is no deeper region around the shear layer, 
where the instability may be amplified by /3-effect, so that we expect the threshold to be much higher and maybe the 
instability to be of a different kind, as the very small azimuthal wave-number observed by [l2| may suggest. 

We reported on figure El the data found in his paper. To check the effect of the geometry, he computed the 
basic-state shear layer in a spherical shell geometry, then took the profile obtained at mid-depth and extended it to 
all z, generating a z-invariant profile in a full sphere [seell2t §5]. He finally computed the instabilities of that flow. By 
artificially removing the inner sphere (and by that way loosing the singular geometrical property described above), 
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FIG. 12: Critical Rossby numbers obtained by iHollerbacbl (2003, summary of his fig. 4 and 8a) showing the impact of the 
slope on the stability threshold. We also plotted the corresponding results obtained with our QG-model (basic state shear layer 
obtained with an inner core, instabilities computed without) including the critical mode. When removing the inner core (ic), 
the threshold drops, and positive and negative Rossby numbers get closer. Our QG-model is close to the fully 3D results of 
Hollerbach 



IHollerbacbl recovers stability thresholds very close to the Ro > case, which is much less modified by this process, 
for it is anyway controlled by the moderate slopes of the outer sphere. 

We repeated exactly the same procedure with our QG-model : wc computed the basic state shear layer with an 
inner-core (ro = 0.35), and then computed the stability threshold of this vel ocity profile in a full sphere geometry 
(without inner core). The results are compared with the 3D computations of IHollerbacbl on figure For Ro < 0, 
the QG-model gives thresholds about 30% higher, and only 6% higher for Ro > 0. The critical modes obtained are 
also in very good agreement. 

This geometrical singularity also prevents any attempt to compute the instabilities of the spherical shell setup with 
a QG-model. We used the dynamo benchmark code described by [(| to perform a fully 3D numerical calculation at 
E = 10~ 4 ' 5 . The critical Rossby number we find Ro c = 0.16 is in very good agreement with the experimental one 
and with the numerical results of [f3 | as shown in figure ^ 



C. Experimental study 

1. Experimental setup 

Figure shows a sketch of the experimental set-up we use to study the onset of instability of the Stcwartson layer. 
It consists of a plexiglass ellipsoid put in rotation by a brushless motor along its minor axis with a controlled angular 
velocity Qq. The working fluid is water. On the axis of rotation a shaft associated to a step motor can drive cither 
an inner core or two disks. Exact dimensions are reported in fig. 1131 and further details can be found in |l9j . 

To determine the stability threshold, we use a flake solution (Kalliroscopc AQ 1000) and illuminate our container 
with a light sheet in a plane including the rotation axis. Variations of brightness allow us to visualize the shear-field in 
the illuminated plane. When the shear-field is steady, we assume that the flow is stable, and when the shear becomes 
unsteady (oscillations in time), the flow is unstable. Two pictures of the observed patterns are shown on figure ITU 
showing that the features of the flow are aligned with the rotation axis. The thin bright and dark columns appear 
at the onset and are blinking, showing the advection of the pattern of alternating vortices. Our visual determination 
shown in figurc^is in good agreement with the numerical calculations of |l2j . This makes us confident in the reliability 
of the optical method. 



2. Experimental results 



The threshold is obtained for different values of the Ekman number with the two geometries (inner core or two disks) 
for prograde differential rotation. Results for the stability threshold are summarized in figure 1151 and quantitative 




FIG. 13: Sketch of the experimental setup with disks geometry. Spherical shell geometry is obtained by replacing the disks by 
a sphere (4cm radius) on the inner rotating axis. 




FIG. 14: Instabilities near the threshold as seen in our experiments (for E ~ 10 °). The vertical shaft is along the rotation 
axis. The features shown by the flakes are vertical, supporting the QG-approximation for our model. 

data are reported in tables IC 31 and IC 41 

Note that in the case of negative differential rotation, we were not able to determine accurately a stability threshold. 
In some experiments hysteresis appeared, correlated with the apparition of air bubbles in the tank. We have also 
observed very low frequency oscillation is some cases. We do not completely understand these results. For the disks 
the data are not even on a straight line. This may be correct, as the numerical simulations also show some oscillations, 
due to the integer nature of the azimuthal wave number (see fig. I15|) . 

The disk case seems to agree with the E 1 ^ 2 scaling, whereas the sphere case is significantly steeper. However with 
only one decade of data, and for such high Ekman numbers, we do not expect the asymptotic laws to hold. 

The numerical calculations for the disk model (see also i jV C 2|) . also shown in figurc^]try to mimic the experimental 
setup : it is a split-sphere with flat polar caps with no height jump (we cannot model discontinuities in the height of 
the container). 

The spherical-shell experiments are in agreement with the numerical calculations of |l2j . as shown on figure ^ The 
strong difference between positive and negative differential rotation is explained by the geometrical properties of this 
setup, as discussed above. 
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FIG. 15: Experimental determination of the stability threshold for two geometries, compared with the numerical calculations 
using a quasi-geostrophic model, sphere 3D is computed with a fully 3D model. The negative differential rotation data are not 
well defined. The numbers are the slopes of the linear fits. 

VII. CONCLUSION AND DISCUSSION 

In this study, we give some insight into the instability of a Stewartson shear layer in the very low Ekman regime. The 
QG-model with the pseudo-stream function formalism is designed for finite slopes as it handles the mass conservation 
correctly. Despite failing to describe the E 1 / 3 -layer, it has proved to be a valuable tool for numerical experiments 
to study very low Ekman numbers. It can also be used for non-linear calculations as long as the quasi-geostrophic 
hypothesis stands, and may be thus useful for other ocean and planetary interior models. 

At very low Ekman numbers E < 10 -6 , we report an interesting radial structure, comparable to the one exhibited 
in the thermal instability case by @- However, if the global structure of the flow is quite similar to the spiralling 
thermal convection flow in a rapidly rotating spherical shell, the radial extension of the instability is surprisingly 
independent of the Ekman number. 

Our model allows us to compare flat top and bottom geometries with varying depth containers. Two different 
stabilizing processes are at work, leading to different balances and scalings for the threshold : for flat containers the 
Q E 3 / 4 scaling applies, whereas with a slope a (3E 1 / 2 scaling is predicted. This is supported by numerical simulations 
for different geometries, and the following table summarizes the different cases for the instability of a Stewartson layer 



geometry Ro c sign of Ro 



f c comment 



j3 = E 3 ' 4 Rot = -Ro~ E- 1 / 4 ? Re c = Cte, Kelvin-Helmholtz-like 

constant (3 ^ fiE 1 / 2 Rof ~ —Ro~ E^ 1 ^ E 1 ^ curvature effect, Rossby wave 

variable (3 [3E 1 / 2 Ro+ ^ — Ro~ E~ X I A E 1 ^ spiralized Rossby wave 

height jump ? not symmetric E~ x l 4 ? QG-model breaks down 

|T^ studied the difference between positive and negative Rossby numbers, and pointed out that the geometry is the 
crucial feature, through the radial derivative of the height of the container. The numerical calculations he made shed 
some light on the dramatic changes in the critical mode when switching from positive to negative Rossby numbers 
observed in some real experiments. However, ignoring the Rossby wave processes at work in varying depth containers, 
he misunderstood his numerical results. 

As we demonstrated in this paper, all the geometrical effects are due to a fundamental symmetry breaking in the 
Rossby wave propagation mechanism : depending on the sign of Ro and /3, the instability will develop on one side of 
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the shear layer or the other. This symmetry breaking happens with any geometry, without the need for abrupt changes 
in the height of the container (even though this asymmetry will be magnified by such discontinuities). Furthermore, 
the spherical shell setup is singular : when Ro < 0, the Rossby wave instability cannot develop on either side of the 
shear layer, leading to the singular results observed in this configuration. Within this framework, the recent results 
of 01 ma y be interpreted correctly. 

What still remains unclear is the nature of the instability in that particular case (spherical-shell and Ro < 0). Is 
it just a Rossby wave that needs much stronger forcing to grow, or another type of instability with significant 3D 
effects ? This is even more annoying because this geometry is of geophysical interest for gaseous planets or planetary 
interiors. It is also the geometry used in the sodium experiment presented by @, aimed at the study of induction 
processes in this type of flow, which are suspected to be able to drive a dynamo. In addition, the non-linear and 
turbulent regimes have to be investigated, by means of numerical calculations and quantitative velocity measurement 
in experiments. It is clear that we need experiments at lower Ekman numbers, especially for spherical shells. 
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APPENDIX A: EKMAN PUMPING IN A SPHERE 

The Ekman pumping flow is given by [iH ] 



E 1 ' 2 „ (nAu±u\ , A1 . 



n.e z 

with n the outward unit vector normal to the surface and u the velocity jump between the geostrophic interior and 
the boundaries. This Ekman pumping flow is the next order term in the E 1 / 2 development, and is deduced from the 
quasi-geostrophic, 0-order flow. However, it has a dynamic action on that flow, and is thus important for the dynamics 
of the system. For simplicity we assume that the boundaries are at rest in the rotating frame. It is straightforward 
to apply the result to moving boundaries by replacing the velocity by the jump in velocity. 



1. Spherical coordinate system 



We can separate this 



l E l/2 

2 



n.V A 



n A u 



±n. V A 



\ 



J 



We now use the spherical coordinate system (p, 9, <j>) to develop that expression, because it is the natural one for that 
boundary layer (n = e p ). Term a reduces to 

_i/ 2 / |tanfl| _ du p _ 2 



and term b can be written as 



a=|cos#| ' ± — ■ Ug u 

2p dp p 



6=|cos6'| V 2 l j- 1-""" i u ^ _|_ ^ 



tanf 



On the boundary, we have p = 1, and so 



E 1 ' 2 ( Itan 



u.n\ ±L = 



2 |cos6»| 



^l 1 / 2 V 2 



2p 



(±u g - u^) ±uj p - 2u p - 



(A2) 
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2. QG-model flow 

We now assume the use of our QG approximation described in ||n] We need to translate IA2I into cylindrical 
coordinate system (r, </>, z). We have 

p sin 9 — r p cos 9 = z 
ug = u r cos 9 ~ u, sin 9 u p = u r sin 9 + u z cos 9 



d dr d dz d , . d n d 

TT- = a o - + "o o = SU10— + cos6l -5- 
op op or op az or az 

Furthermore, u r and are z-independent, and u z is a linear function of z (see section Til CJI 



We find 



We also need to develop lo p 



where cu r is given by 



and finally 



— v 

u z = z ?u r (r, 6) (A3) 

1 — r z 



du„ -2r 

; u r (A4) 



dp 1 — r 2 r 



u! p = lo z cos 9 + u> a sin ( 



1 du z — z du r 
r r 86 1 — r 2 cM 



£i/ 2 ( -r 5 r 



l±£ 2(l-r 2 )V4 ^2(1-^)1/2^+2 1-^^ 



This expression can be treated just like the 0-order no-slip boundary condition : 



u.n \, T = ru r 



]±L — i a r ± y/l — r 2 u z 
from which we deduce u z at the top and bottom boundary 

±u z \ ±L = (1 -r 2 ) (u.n| ±i -ru r ) 

and from the linear z-dependence of u z we obtain 

du z E 1 ' 2 ( r (du r 1 \ 5r 

~dz~ ~ 2(1 -r 2 ) 3 / 4 + I - T 2 ^00 ~ 2 U0 J ~ 2(1 -r 2 f/ 2 '' 

— r 



' 1 - r 2 ' 

E 1/2 .P(u r , Ul p,r) + (3u r 



(A6) 



It has to be emphasized that this results from a development in E 1 ! 2 , which is no more valid when approaching the 
equator. 
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3. Validity 

In the expression I A7I we have an order E 1 / 2 and order 1 term. The development is no more valid when these two 
terms are of the same order of magnitude. Setting r = 1 — e, we have 1 — r 2 ~ 2e. Keeping the dominant terms in 
1/e, the limit of validity is given by 

E 1 ' 2 5 1 1 

U r 



2(2e) 3 / 4 2 (2e) 3 / 2 2e 
from which we get 

i fsV'" 



£ 2 / 5 (A7) 



which shows the i? 2 / 5 singularity on the equator [seeOOJ. The angular extension is easily obtained from geometrical 
considerations : 

s = cos 59 ~ 1 - 56 2 /2 ~ 1 -e 

so that 56* = (2c) 1 / 2 and 

<56» - S 1 / 5 (A8) 



APPENDIX B: A SIMPLE MODEL FOR THE RADIAL SCALING 



1. Approximations and assumptions 

In figure 0abc the perturbation of the shear layer spreads across the whole volume whereas in the flat container 
case (fig. or even the thermal convection case in a sphere [see@, it is radially localized. In this section we derive 
a simple heuristic model that predict all these properties outside the Stewartson shear layer. 

To simplify our study, we make a few approximations : 

• we neglect the curvature; 

• we drop the viscous dissipation term, but keep the Ekman friction; 

• we assume small forcing, ignoring the local vorticity term. 

The curvature may be negligible if the scales involved are small enough. Looking at the results from section IV Al we 
expect the scales to shrink while lowering the Ekman number. Thus, we replace the angular variable <f> by y. Dropping 
the viscous term lowers the order of the differential equation and makes it easier to handle. Numerical simulations 
with no viscous term or an artificially weakened viscous term show that the structure stays very similar. In addition, 
the threshold is expected to be only slightly affected in the asymptotic regime. This is because the Ekman friction 
is also a dissipative term. Finally, we can neglect the local vorticity term because in the linear stability case it is a 
higher order effect, as shown in section fVBI We end up with 

„ / tt 9cj <90\ „ / du z \ 

m +Ro { U dy- + u ^)= 2 {-df) ™ 

The expression of the Ekman friction involves several terms. We will only keep the vorticity contribution, so that 

H^S-pUr-aE 1 ^ (B2) 

where a is the geometrical factor for the Ekman friction term, and for a sphere a = 1/2(1 — r 2 )~ 3 / 4 (see eq. IA6I) . 
We now assume a perturbation of the form 



*(r,y,f) - i>{r)e i{ - kv - ut) 
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and wc obtain 



(-iu + RoUik + E^ 2 aj (k 2 - V + {^°~^~ ~ 2 (^j = ( B3 ) 



In a more convenient form, we can write 

1 C2/9- Rn§B\(t,> - RnTIk-iK 1 / 2 ^ 

ip (B4) 



d 2 ip 2 / 1 (2/3 - itof£) (w - RoUk - iE^a) 
<9r 2 _ I k fa- RoUkf + Ea 2 



2. Outer part, large (3 



In the outer part the basic state velocity profile U is vanishing, as well as ^p. This simplifies the previously 



established relation IB4 



dr 2 



k 2 1 + =f- 



2(3 lu - iE 1 ' 2 a 



k u) 2 + Ea 2 



(B5) 



In the case of small Ekman number, Ea 2 is negligible compared to to 2 ~ E 1 ' 2 , but we have to keep the small imaginary 
part. We now set 



If 1 2(3 1 3> 1, we can rewrite IB5I as 



7 = k 



d 2 ^ 2( 



2/3 
few 



1/2 



2i/3£ 1 / 2 a 
/few 2 + 2(3uj 







(B6) 



and the solution for constant a and (3 is of the form 



?/> = cxp ±7- — — r exp(±i7rj 



ku, 2 + 2(3uj 

This implies a rapid oscillation of typical length scale 7 

7^-^31 

modulated by a slow exponential decay with typical length scale S given by 

S ~ 7~ 1 iT 1/2 a- 1 w 



(B7) 



(B8) 



Inside this envelope, the radial length scale 7 of the perturbation decreases when (3 increases. It corresponds to the 
spiralling of the Rossby wave observed for variable (3 (see fig. \^jpc) 

Using the scalings for the Stewartson layer instability k ~ i?" 1 ' 4 and lo ~ E 1 / 4 , relation IB8I leads to 



*~(2| 



)-V3 a -l£0 



(B9) 



which implies that the radial extension of the perturbation is independent of E. This is a fundamental difference with 
the thermal convection case investigated by [24| and 0] where the scaling k ~ _E -1 / 3 and uj ~ S 1 / 3 implies 



J - (2| / 3|)- 1/2 a- 1 E 1 / 6 
which corresponds to a slowly decaying radial extension when decreasing E. 



(BIO) 
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3. Inner part, small f3 



The inner part behaves exactly as the outer one, replacing lu by uj — RoUk. For a spherical shell it corresponds to 
the small j3 case. 

If 1 2(3 1 < kw, we can rewrite lB5l as 

d 2 ip 2 2i(3E 1 / 2 a 



w-^V-T^Tw)* (BU) 

which actually exchanges the length scale of the exponential decay and of the oscillation compared to the outer region. 
The solution for constant a and f3 is then 

/ j3E 1 / 2 a 
^ = exp(± 7 r) exp^7 fca;2+2 ^ r 

which is an exponential decay of characteristic length 7 -1 , inside which there are much larger length scale oscillations, 
which cannot be detected. 

For the Rossby wave, we have u> ~ fc _1 , so that 7 ~ k, with a correction independent of E. This shows that in the 
small j3 case the critical mode has a radial extension comparable to its azimuthal scale. 



APPENDIX C: STABILITY THRESHOLD DATA 

1. Split-sphere stability threshold 

Linear stability threshold values for different Ekman numbers E. NR denotes the number of radial grid points used 
for the calculation. These values where obtained by QG numerical calculations with no-slip boundary conditions and 
including Ekman friction. 
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2. Boundary conditions and viscosity dependence 

Comparison between Ro c of different numerical models for the split-sphere geometry. Values are the Ro c obtained 
by a model divided by the reference Ro c (ns,ep) for Ro > at a given Ekman number. The models are referenced 
by ns and fs respectively for no-slip and free-slip boundary conditions, ep for the Ekman pumping, dv for horizontal 
divergence and vs for bulk viscosity; x2 stands for calculation with doubled number of radial grid points. A denotes 
the absence of this feature. 
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3. Spherical shell experimental threshold 

Stability threshold experimental data obtained for the spherical shell geometry, in water (y = 1CP 6 m 2 /s). See 
aVICll for details about the experimental setup. 



E 




stable Ro 


unstable Ro 


3.51 10- 


'S 


0.175 


0.183 


2.30 10" 


-5 


0.135 


0.141 


1.21 10" 


■5 


9.29 10~ 2 


9.56 10~ 2 


7.30 10" 


■6 


6.93 10- 2 


7.10 10~ 2 


4.88 10" 


-6 


5.30 10~ 2 


5.52 10~ 2 


3.67 10" 


-6 


4.48 10" 2 


4.56 10~ 2 


2.30 10" 


-5 


-0.302 


-0.313 


1.21 10" 


-5 


-0.232 


-0.246 


6.09 ID- 


■6 


-0.165 


-0.179 



4. Disks experimental threshold 

Stability threshold experimental data obtained for the disks geometry, in water (y = 10~ 6 m 2 /s). See WI C II for 
details about the experimental setup. 
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